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Abstract 


We are concerned with the classification of positive radial solutions for the system Am = m^. 
Am = /(|Vm|), where p > 0 and / S C^[0,oo) is a nondecreasing function such that f{t) > 0 
for all t > 0. We show that in the case where the system is posed in the whole space 

-p/(2p+l) 

ds = oo. This is the counterpart 


such solutions exist if and only if 


{J F(t)4ty 


of the Keller-Osserman condition for the case of single semilinear equation. Similar optimal 
conditions are derived in case where the system is posed in a ball of . It , q > 1, 

using dynamical system techniques we are able to describe the behaviour of solutions at infinity 
(in case where the system is posed in the whole or around the boundary (in case of a ball). 


1 Introduction 


In 1950s Keller [19] and Osserman [25| obtained independently optimal conditions for the existence 
of a solution to the boundary blow-up problem 


Am = f{u) in 0, 

u = oo on dn. 


( 1 . 1 ) 


where D C is a bounded domain and f ^ [0, oo) is a nonnegative increasing function. The 

condition on the boundary dQ in (|l.ll) is understood as limj;-).^;^ u{x) = oo for all xq £ dfi. Keller 
and Osserman obtained that (11.11) has C^(fi) solutions if and only if 



< oo 


where F(s) 



( 1 . 2 ) 


Interestingly, condition (|1.2|) also appeared in other circumstances: it is related to the maximum 
principle for nonlinear elliptic inequalities. For instance, if u G C^(fi) is nonnegative and satisfies 
< f{u) in D, then, if u vanishes at a point in D, it must vanish everywhere in D. We refer 
the reader to Vazquez m and to Pucci, Serrin and Zou [13 CD EE] for various extensions of this 
result. 

Problems related to boundary blow-up solutions have a long history and they can be traced 
back to at least a century ago when Bieberbach |3| investigated such solutions for the equation 
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Au = e“ in a planar domain. Since then, many techniques have been devised to deal with such 
solutions (see, e.g. [T3l [HI [29] for an account on the progress on this topic). Boundary blow-up 
solutions for semilinear elliptic equations with nonlinear gradient terms have been only recently 
investigated (see for instance [I1[1[II1[23]). 

In this paper we investigate a semilinear elliptic system featuring a mixture of power type 
nonlinearities and nonlinear gradient terms. More precisely, we shall be concerned with 


Au = in n, 

Av = /(|Vrt|) in 12, 


(1.3) 


where 12 C is either a ball centred at the origin or the whole space, p > 0 is a real number and 
/ £ C^[0, cx)) is a nondecreasing function such that f{t) > 0 for all 2 > 0. Our study will assume 
that u and v are positive radially symmetric solutions of (II.3p . Note that we do not assume a 
priori any condition at the boundary for neither u or v but this will be needed in the course of our 
analysis as we shall be concerned with the classification of all solutions to (11.31) . 

If 12 is a ball, system (II.3j) was first considered by Diaz, Lazzo, and Schmidt in [7|, in the 
case p = 1 and f{t) = t^. Such choice of exponent p and function / is related to the study of 
the dynamics of a viscous, heat-conducting fluid. The authors in [7| obtained the existence of one 
positive solution and, in case of small dimensions, of one sign-changing solution that blows up at 
the boundary. Their study was further extended to time dependent systems in Diaz, Rakotoson, 

and Schmidt mm- 

We shall first be concerned with the case where 12 is a ball. In such a situation we obtain that 
(|1.3p admits positive radially symmetric solutions {u, v) such that u ot v (or both) blow up around 
912 if and only if 



ds 


p/(2p+l) 


< oo. 


(1.4) 


^ Jo 


This can be seen as the analogous condition to (II.2p obtained by Keller |18] and Osserman |24] 
for (|l.ll) . We also provide a complete classification of radially symmetric solutions in such a case. 
Moreover, we shall obtain (see Theorem 12.41 below 1 that the equation 


A'^u = f{\Vu\) in Br 


has (not necessarily positive) radially symmetric solutions that blow up at the boundary OBr if 
and only if 



If /(2) = 2'^, q > 1, we are able to give the exact rate at which the components u and v blow up 
at the boundary. In such a setting we use dynamical systems tools for cooperative systems with 
negative divergence. 

Further, condition (|1.4I) appears again in the study of (|1.3p in the case 12 = Again when / 
is a pure power type nonlinearity we shall be able to precisely describe the behaviour of solutions 
at infinity. 
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2 Main results 


Let us first present the analysis of system (II.3p in the case where is a ball. Namely, we shall first 
investigate the system 

f Au = in Bfi, 

{ 2.1 

[ Au = /(|Vri|) m Br, 

where Br C N > 2 is the open ball of radius R > 0 centred at the origin, p > 0 and 
/ G (^^^[0, oo) is a nondecreasing function such that f{t) > 0 for all t > 0. Let F be the antiderivative 
of / that vanishes at the origin (see (jl.2p L Sometimes in this paper we shall complement the system 
(j2.ip with one of the following boundary conditions: 

• either u and v are bounded in Br; 


or 


or 


u is bounded in Br and lim v(x) = oo; 

\x\^R 


lim u(x) = lim vix) = oo. 
\x\/'R kl/'R 


( 2 . 2 ) 


(2.3) 


From the first equation of (12.11) it is easy to see that the situation \\vd.\^^\^^ru{x) = oo and v is 
bounded in Br cannot occur. 


Theorem 2.1. We have: 

(i) All positive radial solutions of \2.1\} are bounded if and only if 

ds 

~ Ts 


i: 


p/{2p+l) 


= OO. 


(2.4) 


(a) There exists a positive radial solution {u,v) of 112.1\) satisfying 112.2\) if and only if 

sds 


i: 


p/{2p+l) 


< OO. 


(2.5) 


(Hi) The exists a positive radial solution {u,v) of i2.1\) satisfying 112.3\) if and only if 

ds , r°° sds 


/ 


{[rmt) 


p/{2p+l) 


< OO and 


p/{2p+l) 


= OO. 


( 2 . 6 ) 


By taking /(t) = e* and estimating the integrals in (|2.4l) - (|2.6p we find: 
Corollary 2.2. Consider 

{ Au = in Br, 

Av = in Br. 

Then any solution of 7| ) is either bounded or satisfies i2.2\) . 

We now let f{t) = t‘^, q > 1. From Theorem 12.II we obtain: 


(2.7) 
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Corollary 2.3. Consider 


in Br 
in Br 


( 2 . 8 ) 


j Au = 

I Av = m Or, 

where p > 0 and q>l. Then we have: 

(i) All positive radial solutions of \2. are bounded if and only if 

, 1 
p <1 and !<(?<-. 

P 

(ii) There exists positive radial solutions of I12.8\) satisfying 112.2^) if and only if 

1 ' 


q>2(l + 


P 


(Hi) There exists positive radial solutions of 112.8\) satisfying 112., if and only if 

1 / 1 

— — 

P ^ P 


The three regions A, B and C in the pg-plane that correspond to the cases (i), (ii) and (iii) in 
Corollary 12.31 are depicted below. 



Figure 1: The three regions described by Corollary 12.31 


In the two pictures below we used MATLAB to plot the solution (tt, v) of system (|2.8I) for g = 3 
and p = 2, p = A and for various space dimensions N. 

Our next result deals with the biharmonic problem that derives from ()2.ip by taking p = 1. In 
this case we are able to deduce optimal conditions for the existence of a boundary blow up solution. 

Theorem 2.4. Let R > 0. The problem 


A^u 


/(|Vu|) 


u = oo 


4 


in Br, 
on OBr, 


(2.9) 














ds 


( 2 . 10 ) 


has radially symmetric solutions if and only if 






Figure 2: The solution u (left) and v (right) for system (I2.8p with boundary condition (j2.2p in the 
case p = 4, g = 3 and for various space dimensions N = 2, N = 20 and N = 40. The graph of the 
u-component is restricted on the vertical axis to the interval [0,50]. 




Figure 3: The solution u (left) and v (right) for system (12.8p with boundary condition (j2.3p in the 
case p = 2, q = 3 and for various space dimensions N = 2, N = 20 and = 40 by restricting the 
vertical axis to the interval [0, 50]. 
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Remark 2.5. (i) In Theorem 12.41 we do not require neither u nor Au to be positive in 

(ii) Classification of radially symmetric solutions for A'^u = may be found in [8l 1201 [211122j . 

We shall next be interested on the behaviour at the boundary of solutions to (j2.8p that satisfy 
either (j2.2|) or Consider 


' Au = v^, u > 0 in Br, 

< Au = iVul"^, u > 0 in Br, (2-11) 

V = oo, on OBr. 

Note first that according to Corollary 12.31 the system (j2.1ip has radially symmetric solutions if and 
only if g > 1/p. 

Our main result regarding the behaviour of the radially symmetric solutions to (j2.11|) is as 
follows. 

Theorem 2.6. Assume p,q > 1, {p,q) ^ (1)1) OLud let {u,v) be a positive radially symmetric 
solution to \2.11\) . Then 


g+2 


lim (R— = 

\x\RR 


(1 +2p)g(g + 2)(g+pg+ 1) 
(pq — 1)2+9 


1 

pq-1 


Also, 


(i) If q > 2(1 + 1/p), then there exists L := G (0,oo) and 


lim 


L-u{x) \{l + 2p){q + 2)'P{q+pq+l)'P 


pq-2(l+p) 


{R-\x\)^^T^ 
(ii) If q = 2(1 + 1/p) then 


{pq - l)2p+i 


1 

pg-1 


lim 

(Hi) If q < 2(1 + 1/p) then 

2 + 2p — pq pq — \ 

lim (R — Ixl) u{x) = - 

\x\/'R' 2 + 2p-pq 


u{x) l{l + 2p){q + 2)P{q + pq + 1)P 


{pq — l)2p+i 


pq-l 


(1 + 2p){q + 2)P{q+pq + 1)^* 
{pq - l)2p+i 

Now we shall be interested in the system (jl.3p posed in the whole , namely 

f Au = yP 


pq-l 


m 


„N 


[ Av = /(|Vtt|) in M^. 

Our main result in this case is as follows. 

Theorem 2.7. ITe have: 

(i) The system 12.16]) has positive radially symmetric solutions if and only if 

ds 


/ 


{[F(m) 


p/{2p+l) 


= OO. 


( 2 . 12 ) 


(2.13) 


(2.14) 


(2.15) 


(2.16) 
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(a) Assume f{t) = , where q>l>p and pq < 1. Let {u,v) be a positive radially symmetric 

solution. If 

P{<1 “ < 2{N-1) (2.17) 


then 


lim 2 ± 2 £^ 

1^1 ^CXD ^ 1—pq 


1—pq 


u{x) _ \{2 + q){N{l - pq) + p{2 + q)y/P{N{l - pq) + q{2p + 1)) 


2p+2 — pq pq~l 

(1 — pq) p (2 + 2 p — pq) p 


P 

pq-1 


and 


lim . „ 

II Q'\~2 

|ai| ^OO ^ 1—pq 


v(x) r(2 + g)(iV(l - pq) + p(2 + q))'^(N(l - pq) + q(2p+ 1)) 


(1 — pq)'^+i 


pq-l 


Remark 2.8. (i) Condition (12.171) hold in particular if 5 < 2. 
(ii) It is easy to see that (uo,vo) given by 


uo(x) 

< 

Vo(x) 


(2 + g)(iV(l _ pq) + p(2 + q))^/P(N(l - pq) + q{2p + 1)) 1 


2p+2—pg pg—1 

{1—pq) p {2 + 2p — pq) p 
{2 + q){N{1 - pq) + p{2 + q))^{N{1 - pq) + q{2p + 1)) 


(1 


pq-l 


2 + 2p — pq 

x\ i-p? 


1-pq 


is a solution of (12.161) with f{t) = that vanishes at the origin. Theorem 12.71 states that any 
positive radial solution {u,v) of (12.161) with f{t) = t^ behaves like (uo,vo) at infinity. 

The remaining of the paper is organised as follows. Section 3 contains a detour in dynamical 
systems; here we state the main tools which we use to study the asymptotic behaviour in Theorems 
12.61 and O The following sections contain the proofs of our main results. 


3 A detour in dynamical systems 

For any points x = (xi, X 2 , X 3 ), y = {yi,y 2 ,y 3 ) in dehne the open ordered interval 

[[x, y]] = {z G : X < z < ?/} C 
Consider the intial value problem 

jCt = g{C) forteM, 

\ C(0) = Co, ^ ^ 

where y : — >■ M is a function. This implies that for any Co £ there exist a unique solution 
C of (j3.1j) dehned in a maximal time interval. We denote by 4>{-,Co) the flow associated to (j3.ip . 
that is, t I—)■ f){t,fo) is the unique solution of (|3.ip dehned in maximal time interval. We shall 
assume that the vector held g is cooperative, that is 

do ■ 

^ > 0 for 1 < i,j < 3, z 7 ^ j. 

The following results are due to Hirsch [MKH]. 

Theorem 3.1. (see |17l Theoreml]) Any compact limit set of 13.1\) contains an equilibrium or is 
a cycle. 
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Definition 3.2. A circuit is a finite sequence of equilibria Ci, ^ 2 ; • • •; Cn = Ci; ^ 2) such that 
n VK^(Ci+i) is non-empty, where W^, denote the stable and unstable manifolds. 

Remark 3.3. If all equilibria are hyperbolic and their stable and unstable manifolds are mutually 
transverse, then there cannot be any circuit. 

Theorem 3.4. (see m Theorem 2]). Let K be a compact set such that: 

(i) All equilibria in K are hyperbolic and there are no circuits. 

(ii) For any T > 0, the number of cycles in K having period less than or equal to T is finite. 
Then: 

(a) Every limit set in K is an equilibrium or cycle. 

(b) The number of cycles in K is finite. 

Theorem 3.5. (see [Ml Theorem 7]) Let Ci, C 2 £ 1^^ such that C,i < C 2 - If 

divfif < 0 m[[Ci,C2]], 

then as. 1\) has no cycles in [[Ci)C 2 ]]- 

Definition 3.6. A subset A C is said to be positively invariant for the flow (j) if (/>(f, (() £ A for 
all C S A and f > 0 . A is called invariant for f if 

4>{t, A) = A for all t > 0, 

that is, for any z £ A and f > 0 there exists C, £ A such that <i){t,C) = 

The following notion of chain recurrence is due to Conley mM- 
Definition 3.7. Let A C be a nonempty positively invariant subset for cf and C G 

(i) For e > 0 and f > 0 , an (e,f)-chain from (( £ A to C' £ A is a sequence of points in A, 
C = Ci,C 2 , ■ ■. ,Cn,Cn+i = C' and of times ti,t 2 , ■ ■ ■ ,tn > t such that |0(L,Ci) - Ci+i| < £• 

(ii) A point C € A is called chain recurrent if for every e > 0, f > 0 there is an (e, t)-chain from f 
to C in A. 

(iii) The set A is said to be chain recurrent if every point C £ ^ is chain recurrent in A. 

Theorem 3.8. (see [H Theorem 3.3], [24[ Lemma 1.4]) If A is connected and chain recurrent, then 
for any (,(' £ A and any £,t > 0, there exists an (e,f) chain from C, to C,'. 

Consider now the initial value problem 

it = G{t,i) forf> 0 , 

^( 0 ) = 6 e m3, 

where G is a function on (0, 00 ) x M^ such that G{t, •) —>■ 5 uniformly on compact subsets of M^. 
We shall say that (j3.2|) is asymptotically autonomous with the limit problem (j3.1|) . We denote by 
$(•,^ 0 ) the semiflow defined by the initial value problem (13.2D . 

The following result will play a crucial role in the proof of Theorem 12.61 





Theorem 3.9. (see [Ml Theorem 1.8]) Let G and assume the trajectory associated 

to i3.S\) is bounded. Then the oj-limit set a;$(^o) has the following properties: 

(a) w<i>(^o) is nonempty, compact and connected. 

(b) is invariant under the flow of 4> of I13.1\} . that is 

4>{t,u;^{Co)) = w<i,(^o) for all t>0. 

(c) w<i)(^o) attracts that is, 

^o ))0 as t —> oo. 

(d) w#(^o) is chain recurrent for cj). 

In particular, Theorem 13.9f dl states that the invariant set consisting of two equilibria ei, 62 and 
a heteroclinic orbit that joins them, or a homoclinic orbit connecting an equilibrium point 63 with 
itself cannot be the cj-limit set of an asymtotically autonomous semiflow. 


4 Proof of Theorem 12.1 


A useful result in proving Theorem 12.11 is the following lemma. 
Lemma 4.1. We have 

ds r°° ds 




p/{2p+l) 


< 00 if and only if 




2p/{2p+l) 


< 00 . (4-1) 


Moreover, 

poo 

Jl 


ds 


2s \p/(2p+l) 


< 


i: 


ds 


2p/{2p+l) 


< 


i: 


ds 


2s 


p/{2p+l) 


Proof. Let H : [0, 00 ) —> [0, cx)) be given by 
H{s) = 2^ F{t)dt- 


for all s > 0 . 


Then 


Therefore 


H'{s) = 2F{s) — 2^ y \/Y(^dtJ \J f{s) for all s > 0. 

= F{s) - ^/7^ [ y/f(fl}dt = [ f{t)dt - [ \/l{fl)dt 

^ Jo Jo Jo 

- Vf(s)y/f(t)dt = y \/7(i)(\/7M- Vfis))dt 


JO 

< 0 for all s > 0 . 
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Hence, H is nonincreasing which yields H{s) < H{0) = 0. This further implies 

2 J F{t)dt — 'J for all s > 0 . 

ds 


So, 




2j^Fit)dt) 


p/(2p+l) 


> 


i: 


ds 


2p/{2p+l) 


In order to establish the first inequality in our lemma, let h : [0, oo) —>■ M be defined by 

h{s) = J F{t)dt “ ( y V for all s > 0. 

Then, 

^ = F{2s) - y' yW)dt 

> [ f{t)dt - y/ f{s) [ y/J{t)dt 

Js Jo 

= [ fis + t)dt- ^/f{s) [ y/f{^dt. 

Jo Jo 

= j (/(i + s) - 

> 0 for all s > 0. 

It follows that h is nondecreasing which yields h{s) > h{0) = 0 for all s > 0. This implies 


Therefore 


2s 


i: 


p/{2p+l) 


> 


2p/{2p+l) 


for all s > 0 . 


(4.3) 


ds 


2s 


p/{2p+l) 


< 


ds 


2p/{2p+l) ■ 


This concludes the proof. 


□ 


Proof of Theorem 12.11 It is enough to prove (ii) and (hi). We shall divide our proof into three 
steps. 

Step 1: Let {u,v) be any positive radial solution of 112.1\) . Then, letting w = u' we have 


vP{r) 

N 


< w'{r) < v^{r) for all 0 < r < R, 


and 


f{w{r)) 


N 


< v"{r) < f{w{r)) for all 0 < r < ii. 


(4.4) 

(4.5) 


Indeed, we note first that {w, v) satisfies 


w'{r) H- w{r) = v^{r) 

r 


for all 0 < r < ii. 


^-1 ./ 


V ir) H- v'ir) = f(\w(r)\) for all 0 < r < ii, 

r 

l(;(0) = ^'(O) = 0,u(0) > 0. 


(4.6) 


10 



























(4.7) 


< r < R. 


Integrating in the first equation of (|4.6I) we find 

w{r) = f for all 0 

Jo 

This implies that u; > 0 in (0, R) so u is increasing. We now integrate in the second equation of 
(14.61) to deduce 

v'{r) = f t^~^f{w{t))dt for all 0 < r < i?. (4.8) 

Jo 

This means v' > 0 and v is increasing on (0, ii). Using this fact in (|4.7I) we find 

w{r) < —v^{r) for all 0 < r < 72. (4.9) 

Combining (j4.9|) with the first equation of (14.6p we have 

- 1 

v^{r) < w'{r) H- — —u^(r) for all 0 < r < 72 

which implies ()4.4p . In particular re' > 0 in (0,72) so w is increasing. Using (14.8p we obtain 

v'{r) < -^fiw^r)) for all 0 < r < 72. (4-10) 

Using (j4.I0p in the second equation of (j4.6p and also the fact that u) > 0 we deduce the estimate 


Step 2: System 112. 71) admits a positive radial solution {u,v) such that \mir/^Rv{r) = oo if and only 

if 

ds 


i: 


p/{2p+l) 


< OO. 


(4.11) 


Assume first that {u,v) is a positive radial solution of (12.ip with \\m.r/-Rv{r) = oo. 
Using ()4.5I) we have 

v"{r) < f{w{r)) for all 0 < r < 72. 

Multiplying the above inequality by v'{r) and then integrating over [0, r] we have 

v'{t)f{w{t))dt 

< f{w{r)) [ v'{t)dt 

Jo 

< f{w{r))v{r) for all 0 < r < 72, 

which gives 


v'{r){v{r)) < Cf {w{r)) for all 0 < r < 72, where C > 0. 

Multiplying the above inequality by w'{r) and then using (j4.4j) we have 

v'{r )-—^ w'{r)v'{r){v{r))~^/‘^ < Cw'{r)y/f {w{r)) 

for all 0 < r < 72. This further implies 

1 / vP+^/'J‘{r) 

N [ p + 1/2 


< Cw'{r)^f{w{r)) for all 0 < r < 72. 
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Integrating over [0, r] we obtain 

pw[r) 

^P+l/2(^) - t;P+l/2(o) <C ylW)dt. 

Jw(0)=0 


Cw(r) 

lw{0)=0 

Since lim^-^/j ti(r) = oo, we can find p G (0, ii) such that 

\ (2p+l)/2p /■"'W 


/ X (2p+l)/2p , - 

~ ^ J V f{t)dt for all p 


<r<R. 


Using (14.41) we obtain 


w'{r) 


w{r) 


2p/(2p+l) 


< C for all /9 < r < i?. 


Integrating the above inequality over [p, r] we have 

w’{t) 


pr 

J P 


w{t) 


0 


2p/(2p+l) 


dt < C{r — p) < Cr for aW. p < r < R. 


By changing the variable and then letting r R one obtains 

ds 


poo 

Jw{p) 


Hence, 


i: 


2p/(2p+l) 


ds 


< C{R — p) < oo. 


(4.12) 


2p/(2p+l) 


< OO. 


Using Lemma l4.ll this is equivalent to 

poo 

Jl 


ds 


p/(2p+l) 


< OO. 


We now assume that / fulfills (j4.11l) and prove that (12.11) has a positive radial solution (u, v) 
satisfying limr^^u(r) = oo. Looking for radially symmetric solutions of (12.ip we are led to solve 


N — 1 

u"{r) H - u'ir) = v'^ir), r > 0 

r 

TV — 1 

v"{r) H - v'{r) = /(|ii'(r)|), r > 0 


(4.13) 


^ u{r) > 0, v{r) > 0 for r > 0. 

In order to obtain the local existence of a solution, it is more convenient to introduce w = u'. Thus, 
the system (I4.13P reads 

u' {r) = w{r) 


TV - 1 

w'{r) H- w{r) = u^(r) 

r 

v"{r) + ^ ^ ^ v'{r) = /(|u;(r)|) 

w{0) = u'po) = 0,u(0) > 0. 


(4.14) 
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By twice integration, (I4.14p is equivalent to 


/■r 

u{r) = u(0) + / w{t)dt, r > 0, 

Jo 

w{r) = f t^-^vP{t)dt, r > 0, 

Jo 

v{r) = v{0) + f f{\w{s)\)dsdt, r > 0, 

Jo Jo 
^ tt(o) > 0, > 0. 


(4.15) 


Since / is a C^-function, by a standard contraction mapping principle one obtains the existence of 
a solution {u, v) of (14.131) defined in a maximal interval [0, Rmax)- By Step 1, w and v satisfy (14.4p 
and ()4.5I) . Thus, we have 


f{w{r)) < Nv"{r) for all 0 < r < Rmax, 

w'{r) < v^{r) for all 0 < r < Rmax- 


Multiplying the two inequalities in ()4.16p and then integrating over [0, r] we deduce 

F{w{r)) < Nv^{r)v'{r) for all 0 < r < Rmax- 
Multiplying the above inequality by 'w'{r) and using (j4.4p one obtains 

w'{r)F{w{r)) < N{r)v'[r) for all 0 < r < Rmax- 
Fix p € (OjRmax) and denote 

rr 

G{r) := / F{t)dt for all p < r < Rmax- 

Jp 

Integrating (j4.17p over [p, r] and using (14.4p we have 

pw{r) pr 

G{w{r)) = / F{t)dt <N {t)v'{t)dt 

J p J p 

< CK(r)](2p+i)/P 

< G[w'{r)]^'^^^^'>/^ for all p < r < Rmax- 


(4.16) 


(4.17) 


Hence 


G < 


w'{r) 


{G{w{r))^ 
A further integration over [p, r] yields 


pr 

G{r - p) < / 

J O 


p/(2p+l) 


W'{t) 


for all p < r < iZr? 


(G(!u(t))) 


V(5TI) < 


Rn 


By changing the variable of integration and then letting r ^ Rmax we have 

ds 


G{Rmax - p) < [ 
J Vu 


w{p) (n( 


p/(2p+i) ■ 


(4.18) 














Therefore, 


Rm.n.T. ^ C 


i: 


ds 


p/{2p+l) 


< oo. 


We have obtained a positive radial solution (u, v) of (12.11) in satisfying v{r) 

Now, if i? > 0 is any positive radius, we let 


= oo. 


ht) = 


for all t > 0. 


Clearly / satifies (|4.1ip . By the above arguments there exists {u,v) such that 


Au = vP 
Av = f{\Vu\) 


where Br^^^ is a maximum ball of existence. Let 



in Br, 
in Br. 


By taking A = R/Rmax, we deduce that {u,v) satisfies (12.11) in Br. 


Step 3: Proof of (ii) and (in). 

Assume (j2.ip admits a positive radial solution {u,v) in Br that satishes (j2.2p (resp. (l2.3p L By 
Step 2 above, / must satisfy (|4.1ip . Prom (14.3p we have 


(f 


2p/{2p+l) 


for all s > 0. 


Using this fact and working in the same way as we did for estimating (|4.12l) and (I4.18p . there exists 
p € (0, i?) such that 


/•OO 

Jw{r) 


ds 


2s 


< 


F{t)dt^ 


P/(2P+1) Jnj{r) 


f 

J wi 


ds 


2p/{2p+l) 


< Ci{R-r) 


and 


/•OO 

Jw{r) 


ds 


[j Fm) 


p/{2p+l) - r) 


(4.19) 


(4.20) 


for all p < r < i?, where Ci, (72 > 0 are constants. 
Let r : (0, oo) —(0, oo) be defined as 


ds 

r{t) = / -- .. ' for all t > 0. 

Note that P is decreasing and by (j4.11l) we have limi_j.oo r(t) = 0. Prom (I4.19P and p4.20p we deduce 
P(2u;(r)) < C\{R — r) and P(tc(r)) > C 2 {R — r) for all p < r < ii. 
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Since F is decreasing, the above estimates yield 

J 2w{r) > T~^{Ci{R — r)) for all p < r < i?, 

[ w{r) < T~^{C 2 {R — r)) for all p < r < ii. 


Let us recall that 


rr 

u[r) = u{p) + / w{t)dt for all p 

J o 


< r < R. 


(4.21) 


(4.22) 


From ()4.22p and (j4.2ip we find u{r) = oo if and only if 


rR 

/ w{t)dt 
J o 


= oo 


if and only if 


fR 

/ T-^{C{R-t))dt 

Jp 


= oo, 


for some constant C > 0. Hence = oo if and only if 

rC{R-p) 


/ 


F '^{a)da = oo, 


if and only if 


[ T-\a)d, 

Jo 


.a = oo. 


With the change of variable t = T ^(a) we now obtain \\m.r/'Ru{r) = oo if and only if 


i: 


( fo F{t)dt^ 


p/{2p+l) 


ds = oo. 


This implies (12.6|) . 

We now assume that (j2.6l) holds and show that system (j2.ip has a positive radial solution (u, v) 
that satisfies (1231). We proceed as in Step 2. First we obtain the (local) existence of such a solution 
in a ball Br^^^ and then, by the same scaling argument indicated at the end of Step 2 we are able 
to conclude the existence of the desired solution to ()2.1h in Br that satisfies ()2.,‘lh . 


5 Proof of Theorem 12.4 


Let V = Au. Then {u, v) satisfies 


TV — 1 

u’'{r) H- u'{r) = v{r) 

r 

v"{r) + — - -v'{r) = f{\u'{r)\) 

r 

lim u{r) = oo. 
r /'R 


for all 0 < r < i?, 
for all 0 < r < i?, 


(5.1) 


Integrating twice in the second equation of ()5.1h we find that r \ — > ^v'{r) is increasing on 

(0, i?), which further implies v' is increasing on (0, R). Thus, there exist 


L := lim vir) G M U {oo}. 

r/'R 
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If L < oo, then v < L for all 0 < r < ii. This fact combined with a twice integration in (I5.1|) yields 

Lr^ 

u{r) — m(0 ) < for all 0 < r < i?. 

This implies that u is bounded over [0, R) which is a contradiction. Therefore we must have 
\\mr/'Rv{r) = oo. In view of this fact there exists > 0 such that 

v{r) > 0 for all Rq < r < R. 

Let u' = w. Then, from (15.ID we deduce 


N - 1 

w'{r) H- w{r) = v{r) 


for all 0 < r < ii, 


N — 1 

v"{r) H- v'(r) = f(\w(r)\) for all 0 < r < ii. 

r 

From the Hrst equation of (15.21) we find 

= r^~^v > 0 for all Rq < r < R, 

which in particular implies that r i— y r^~^w is increasing on [iio,ii) and so, there exists 

Lq := lim w(r) = lim u'(r) G M U {oo}. 

r /^R r /^R 


(5.2) 


(5.3) 


With a similar argument as above, if Lq is finite we derive that u is bounded which contradicts 
\m\r/-Ru{r) = oo. Hence \\.m.r/-Rw{r) = oo. Integrating (15.3h over [iio,r] we obtain 

v[r) 

L vyijui Oi vyr) / t ui = 

^ Rq 

for all Rq < r < R. This yields 


r^-^u;(r) - R^-^w{Ro) = [ t^-^v{t)dt < v{r) [ t^-^dt = ^ - R(^) , 

Jro JRo ^ ^ ’ 


w{r) < 


v{r)r 


+ 


C 


,N-1 


for all Rq < r < R, 


(5.4) 


where C = C{Ro,N) > 0 is a constant. Since v{r) —>■ oo as r R, we may choose Ri G {Ro,R) 
such that 


C 


< 


v{r)r 


rN-1 - 2N{N - 1 ) 


for all Ri < r < R. 


Combining (j5.4p and (15.5p we obtain 


2N -1 

w(r) < - -v(r)r for all Ri < r < R. 

^ ’ ~ 2N{N - 1) 

We now use this last estimate in the hrst equation of (15.21) to deduce 

v{r) 


(5.5) 


(5.6) 


2N 


< w'{r) < v{r) for all Ri < r < R. 


The same approach is now applied to the second equation of (15.2p in order to deduce 

< v"{r) < f{w{r)) for all Ri < r < R. 


fiw{r)) ^ 


2N 

From now on, we follow line by line the proof of Theorem 12.11 with p = 1 to reach the required 
conclusion. 
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6 Proof of Theorem 12.6 


6.1 More properties of solutions to system (12.lip 

Let {u,v) be a positive radially symmetric solution of (|2.11l) in Br. 
Letting u'{r) = w{r) and v'{r) = 'tp{r) we have 




w'ir) H- wir) = u^(r) 

r 

v'{r) = V'(r) 


■0 (r) H- 

r 

w{0) = 0, u(0) 


^p{r) = 

= m > 0 , ' 0 ( 0 ) = 0 . 


for 0 < r < ii, 
for 0 < r < ii, 
for 0 < r < 12 , 


( 6 . 1 ) 


The next result is a comparison principle between sub and supersolutions of (16.ip which is true 
in virtue of the quasimonotone character of our system. 

Lemma 6.1. Let {vi{r),wi{r),'ipi{r)) and {v 2 {r),W 2 {r),'ijj 2 ir)) he solutions of 


+ -wiir) > Vi{r) 

v'i{r) > 01 (r) 

Q 

0'i(r) + -01 (r) > wl{r) 

^ i(;i(0) = = mi,0i(O) = ui. 


/or 0 < r < 
for0<r< 
/or 0 < r < 


where 6* > 0, 

If 


^2(^) “I— 'W2(r) < Voir) 
r 

V 2 {r) < 02 (r) 

^ 0 

02 (r) + - 02 (r) < wl{r) 

_ 102 ( 0 ) = H2,V2{0) = 7712,02(0) = 122. 
> 0, 77* > 0, ?T7i > 0 , (i = 1, 2). 


/or 0 < r < 12 ^^, 
/or 0 < r < 12 ^“"^, 
/or 0 < r < 


L1>H2 , 7771 > 7772, 1^1 > 1^2, (/7l, 7771,770 0 (^ 2 ,7772, 7^2). 


Then 


( 6 . 2 ) 


(6.3) 


wi{r)>W 2 {r), vi{r) > V 2 {r), 0 i(r) > 02 (r) /or 0 < r <' 7777 ? 7 {i 20 '“^, 12 ™“*}. 

Proof. Without loosing any generality we may assume that ttii > m 2 . Note that the first equation 
in (j6.2l) and (j6.3p can be written as 


^r^u;i(r)^ > r^Vi{r) for 0 < r < 120'“*, 

(^r^W 2 ir)^ < r%P{r) for 0 < r < 120“*. 

Since 01 ( 0 ) = 7771 > 02 ( 0 ) = m 2 , there exist p > 0 such that 

7;i(r) > 02 (r) for all 0 < r < p. 


(6.4) 


(6.5) 
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Using (16.4p and (|6.5p we have 

^r^wilr)^ > r^vf(r) > r^v^ir) > (^r^W 2 {r)^ for all 0 < r < p. 

Hence ^ ^ 

^r^rai(r)^ > (j'^W 2 {r)J for all 0 < r < p. (6.6) 

Integrating (16.6p over [0,r], r < p, we find 

rt;i(r) > 'W 2 {r) for all 0 < r < p. (6.7) 

Using (16.7p in the third equation of (16.2p and (16.3p we get 

^r^'0i(r)^ > ^r®V’ 2 (?’)^ for all 0 < r < p. 

Integrating over [0, r], r < p, we deduce that 

V'i(r) > ' 4 ’ 2 {'I") for all 0 < r < p. 


Let us denote 

R := sup |p > 0 : ui(r) > V 2 {r),wi{r) > W 2 {r),ipi{r) > 'ip 2 i'i") in (0,p)|. (6.8) 

First of all note that i? > 0. Using the same arguments as above we have 


□ 


Lemma 6.2. Let {vi{r),wi{r),'ipi{r)) and {v 2 {r),W 2 {r)^'il) 2 {'f’)) be the solutions of 


iV - 1 

“1- wiir) = v^{r) 

r 

for0<r< 


v[{r) = ifiir) 

for0<r< 


N — 1 

'f’lir) H- - —-01 (r) = wl{r) 

for0<r < 

(6.9) 

t(;i(0) = 0,771 (0) = 7771 > 0,7/71(0) = 0 



— 1 

'W^2(^) 4-7772 (r) = 77f(r) 

r 

for0<r< 72^“^, 


V2{r) = i^2{r) 

for0<r< RJff^^, 


N — 1 

7/7^ (r)+ ^ 7/72(7) = 77;|(r) 

/or 0 < r < 72^“^, 

(6.10) 

V2{Rrn = 00 , 




^ W2{0) = 0, U2(0) = 7712 > 0,'lp2{0) = 0. 
If mi > 7712 then 7?™“^ < . 
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Proof. Let cr > 0 be such that 


pq — 1 

\mi 


< a < 1. 


Since mi > m 2 and pq — 1 > 0, it is always possible to choose such a a. Then 


2p+l g+2 _ pg+g+1 

w{r) = api-^wi{ar), v{r) = api-^vi{ar), ip{r) = a 'il’i{ar) 


( 6 . 11 ) 


satisfies 

( , N -I 

w (r) H- w(r) = v^ir) 

r 

v'{r) = i>{r) 

{ iV-1- (6.12) 

^^(r) H- 'ipir) = w^(r). 

r 

~ g+2 

^ u;(0) = 0, -0(0) = 0, ^(O) = api-^mi > m2. 

Using Lemma [6T] for (^(r),u}(r), V’(r)) and {v 2 {r),W 2 {r),'ip 2 {'>')) it follows that 


v{r) > V 2 {r), w{r) > W 2 {r), > fj 2 {r). 


In particular v{r) = ap‘i-^vi{crr) > V 2 (r). Since ui(cjr) is finite as long as ar < it follows 

^max 

that V 2 {r) is finite as long as r < —. Since is the maximnm time of existence for V 2 {r), it 

follows 


rymax 

^1 ^ j^max 

a ~ ^ 


Using the fact that <7 < 1, we have 


j^max 




^ T^max 

< -Ko , 


which completes the proof. 


□ 


6.2 Proof of Theorem 12.61 

We may always assume R = 1 because if (u, v) is a solution of (12.lip in Bji then 



f 2 + 2 p—pg 

I U{x) = R p-j-i u{Rx) 

1 2+9 

V{x) = Rpi-^ v{Rx) 

for X G i?i, 

for X G Bi, 

(6.13) 

is a solution of (| 2 . 11 |) 

in Bi. In the seqnel we shall assume R = 1. Thus, letting w 

= u' and ij; = v' 

we have that {v, w, if) satisfies 




f N - 1 

w'{r) H- w{r) = u^(r) 

T 

for 0 < r < 1 , 



v'{r) = ' 0 (^) 

for 0 < r < 1 , 


< 

N — 1 

fj'ir) H- ifir) = w^^ir) 

for 0 < r < 1 , 

(6.14) 


u;(0) = 0,u(0) = m> 0,'i/’(0) = 0, 
lim w{r) = lim v{r) = lim 'ip{r) = 00 . 

r r r 
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Further, let 


(6.15) 


«(’’) = ^^( 1 -^)", Kr) = ^{'^-rf, c(r) = ^^( 1 -r)^, 


where 


A = 
B = 
C = 


(1 + 2p){q + 2Y{q+pq + 1 )^ 


{pq — l) 2 p+i 

{l + 2p)l{q + 2){q+pq+l) 


1 

pq-1 


1 

pq-1 


{pq — 1)2+9 
(1 + 2p)i{q + 2)P'^{q+pq + 1 ) 


{pq - l) 9 +P'?+i 


1 

pq-1 


and 


l + 2p q+ 2 q+pq +1 

a = -r, P = -7 = 


pq — I PQ — ^ 

From (|6.14p we deduce that (a, b, c) satisfies 


pq-l 


(1 — r)^a'(r) + ——-a(r)^ = a{V’{r) — o(r)) 
(1 — r)b'{r) + Pb{r) = f5c{r) 

AT 

(1 — r) (^c'{r) H-c(r)^ = '){a^{r) — c{r)) 

a(0) = 0,6(0) = m/B > 0, c(0) = 0. 


for 0 < r < 1 , 
for 0 < r < 1 , 
for 0 < r < 1 , 


(6.16) 


(6.17) 


(6.18) 


We next introduce a new change of variable in the system (I6.18p . by letting r = 1 — e ^ and 
X{t) = a(r), Y{t) = h{r) and Z{t) = c(r) where t = In(j^). Thus, (I6.18h yields 


’ + ^^X{t) = a{YP{t) - X{t)) 

— 1 

Yt = P{Z{t) - Y{t)) 

Zt + ^^Z{t)=^{X‘^{t)-Z{t)) 

_ X(0) = 0, F(0) = m/B > 0, Z{0) = 0 


for 0 < t < oo, 
for 0 < t < oo, 
for 0 < t < oo. 


(6.19) 


The proof of Theorem 12.61 will be divided into three steps. 

Step 1: ^{t) = {X{t),Y{t), Z{t)) is bounded as t —)• oo. 

Let us assume by contradiction that {X{t), Y{t), Z{t)) is not bounded. Then we claim that Y{t) 
is unbounded. If Y{t) is bounded, the first equation of (I6.19P would imply 

Xt + aX{t) < aYP{t) < C for all t > 0. 

which is equivalent to 

(^X(t)e“*y < Ce"* foralH> 0 . 

Integrating the above inequality we easily deduce that X{t) is bounded. Similarly, Z{t) is bounded 
which contradicts our assumption. Therefore Y{t) must be unbounded. 

Let m > m and (u, u), ip) be the solution of (16.11) with the initial conditions 

u)(0) = 0, u(0) = rh, ip{0) = 0 
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defined on the maximum interval {0,R). By Lemma 16.21 we have i? < 1. 
Let {X, Y, Z) be the solution of 


Xt + = a{YP{t) - X{t)) 

— 1 

Y, = f5iZit)-Y{t)) 


Y - 1 


Z{t)=jiXiit)-Z{t)) 


for 0 < t < T := In 
for 0 < t < T, 
for 0 < t < r, 


(t^)' 


( 6 . 20 ) 


, X(0) = 0, y(0) = m/B, Z{0) = 0, 


associated to {v,w,'^). Then {X,Y,Z) blows up at T. 

Since Y is unbounded, we can choose to > 0 such that Y{to) > y(0) = fh/B. Let us set 

X{t)=X{t + to), Y{t) = Y{t + to), Z{t) = Z{t + to). 

Then, one can easily check that {X, Y, Z) satisfies 


Xt + ^^^X{t) > aiYP{t) - Xit)) 

— 1 

for 0 < t < oo 

Yt = l3{Z{t)-Y{t)) 

for 0 < t < oo 

Zt + ilY.zu) > I(x\t) - zm 
— 1 

for 0 < t < oo 


, X(0) > 0, y (0) > m/B, i(0) = 0. 
In virtue of Lemma 16.11 we deduce that 


X{t) > X{t), Y{t) > Y{t), Z{t) > Z{t) 


which contradicts the fact that (X,Y,Z) blows up in finite time. Hence ^(t) = (X(t), Y(t), Z(t)) 
is bounded as t —)■ oo. 


Step 2: Analysis of the autonomous system associated with i6.19\) . 

We shall embed the autonomous system associated to (I6.19P in the whole by considering the 
initial value problem 


Ct = giO for all t e M, 
C(o) = Co G m3, 


( 6 . 21 ) 


where 

( X\ /X\ /a{Y\Yr^-X)\ 

C= y , 5(0=5 ^ = I3{Z-Y) 

\zj \zj \j{x\xr^-z)J 

Using a standard comparison result we have 


Lemma 6.3. Let {X{t),Y{t), Z{t)) be the solution of \6.21\) and to G M. Then 
(i) If X{to), Y{to), Z{to) < 0 then X{t),Y (t), Z{t) < 0, for all t < to- 
(a) If X{to),Y (to), Z{to) > 0 then X{t),Y (t), Z{t) > 0, for all t > to- 


21 















The system (16.211) is cooperative and has negative divergence. It has exactly three eqnilibria, 
namely 0 = ( 0 , 0 , 0 ), 1 = (1, 1 ,1 ) and —1 = (—1, —1, —1). It is easy to check that 0 is asymptotically 
stable. The linearized matrix at 1 and -1 is 


M 


—a ap 0 
0 -/3 /3 , 

IQ 0 “7 


and the eigenvalues Xi, 1 < i < 3, are solutions of 


(A + a)(A + /3)(A + 7 ) — pqafi'y = 0. 
From the definition of a, /3 and 7 in (|6.17l) . we have 


a + l=p/3, /3 + 1 = 7 , ^ + l = qa. 

This shows that Ai = 1 is an eigenvalue of M. Also A 2 + A 3 < 0 and A 2 A 3 = {pq — l)(a/ 37 ) > 0. 
Thus, Re{\ 2 ) < 0, Re{X'i) < 0. So 1, -1 are saddle points with two-dimensional stable manifolds. 
Using Lemma 16.31 and the fact that 0 is asymptotically stable, we deduce that the system (I6.2ip 
has no circuits. By Theorems 13.41 and 13.51 any compact limit set of (I6.2ip reduces to an equilibrium 
point. Hence, any bounded trajectory C) converges both backward and forward in time to one 
of the three equilibria described above. 

Step 3: Analysis of the non-autonomous system 116.19\) . 

Let ^0 = {0,m/B,0) and denote by <k(-,^o) the semiflow associated to (16.191) . By Theorem 13.91 
the w-limit set is invariant under the flow 4> of the autonomous system ()6.2ip . Thus 

(/>(t, w$(^o)) = ‘^'*'(^ 0 ) for all t > 0 . 

Due to the group property of the flow cj), the above equality is true for all t € M because 


= 0(O,w<i.(?o)) = </>(-|t|X|t|,W4,(^o)) = </>(-|t|,W4,(^o))- (6.22) 

Hence, (/>(t,a;$(^o)) = for all t gW. 

Let z G w#(^o)- Since is chain recurrent, for all n > 1 there exist a finite sequence of 

points in cn$(^o) 

^ Cl; C 2 ; ■ ■ ■ ; C/Cn ; Cfcn + 1 ^ 

and a sequence of finite times 

ti, t2, ■ ■ ■, tk^ > n 

such that ^ 

IC>(U, Ci) - Ci+i| < - for all 1 < i < kn- 
n 

In particular, for i = kn we find 

\4>{tkn,CkJ - z\ < - for all n > 1. (6.23) 

n 

Since {Cfc„} C w$(Co) is bounded, it follows that up to a subsequence (still denoted by {Cfc„}) we 
have Cfc„ —>■ Co as n —)■ 00 , for some Co G a;<j>(Co)- By Step 2, (j){t,(o) —^ G {0;1} as t —>■ 00 . 
Using the continuous dependence of the flow (j) on the initial data we can let n —)• 00 in (j6.23p to 
deduce z = i G {0,1}. Thus, a; 4 )(Co) C {0,1}. Since a; 4 )(Co) is connected, by Theorem IS.Ol al it 
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follows that a;$(.^o) = {0} or a;$(^o) = {!}• Assume by contradiction that a;$(.^o) = {0}- Then 
{X{t), Y (t), Z{t)) tends to 0 as t —>■ oo. Then, we may find to > 0 such that 


0 < A(to),T(to),^(to) < 1- 


Let now m > m and {X{t),Y{t),Z{t)) be the solution of ()6.20p . By taking m close to m and using 
the continuous dependence on the initial data of solution to (I6.19p . we may assume 

0 < X{to),Y{to),Z{tQ) < 1 . 

A comparison principle now implies 

0 < X{t),Y{t), Z{t) < 1 for all t > to- 

But the above inequalities contradict the fact that {X{t),Y (t), Z{t)) blows up in finite time. Hence, 
e(t) = (A(t),y(t),z(t))^i. 

Using (I6.15|) it follows, 

u(r)(l — r)^ —> B as r 1, 

tc(r)(l — r)“ A as r ^ 1, (6-24) 

'(/’(r)(l — ry C as r Z' 1. 

Now, the first part of equation (j6.24p implies (I2.13p . Let e > 0. Then there exist 6 G (0,1) such 
that 

(1 — e)A(l — r)““ < w{r) = u'{r) < (1 + e)A(l — r)“" for all r € [5,1). (6.25) 

Assume q > 2(1 + 1/p). Thus a G (0,1). By Corollary 12.31 u is bounded and increasing on (0,1). 

Thus there exists L = lim^y-ii u{r) and integrating (I6.25P over [r, 1], we find 

^ -(!-£)< <T^(l+g) for all rG [5,1). 


1 —q;' ^ (1 — r)^ " 1 — a' 

This proves part (i) of Theorem 12.61 

Assume now q = 2(l + l/p). Thus a = 1. Integrating (|6.25l) over [5, r], where <5 < r < 1, we find 


A (1 -e)< liminf < limsup < A (1 + e). 


Letting e —>■ 0, we get 


In(^) ,./>! In(j^) 

u{r) 


lim ——- 
''/'i ln(^^ 


= A. 


This proves part (ii) of Theorem 12.61 

Assuming q < 2(1 + 1/p), in a similar way as before we derive the proof of part (iii) in Theorem 


7 Proof of the Theorem 12.7 


As in the proof of Step 1 in Theorem 12.11 we obtain that u', v', u, v are increasing and 

- pr 

u'{r) = / s^~^v^{s)ds for all r > 0 , 


v\r) = Z ^ f ^{u^y{s)ds for all r > 0 . 

> Jo 
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This yields 


for all r > 0 


(7.1) 


rv^lO) ,, . rvP(r) 
< u [r) < 


and 


N 

vpg(^0)r^+q 


< v'{r) < 


N 


for all r > 0. 


N<}{N + q) - ' ^ - N 

From (EH) and (I7.2p we deduce that u'{r), v'{r), u{r), v{r) tend to infinty as r —>■ oo. 

Inspired by the change of variables introduced in [18] (see also mm) we define 

,,, , ru'ir) , rv'ir) , rvP(r) , , ru'^ir) 

u[r) v(r) w{r) v'[r) 

where t = ln(r) for r € (0, cx)). A direct calculation shows that {X{t),Y{t), Z{t),W{t)) satisfies 

( Xt = X{Z -{N-2)-X) for all t G M, 

Yt = Y{W -{N-2)-Y) for all t G M, 

Zt = Z{N +pY -Z) for all t G M, 

Wt = W{qZ -qN + q + N -W) for alH G 


(7.2) 


(7.3) 


By L’Hopital’s rule we have limi_j.oo Ai(t) = 2 — + limf_j.oo Z{t). Thus, it is enough to consider 

the last three equations of (17.31) . namely 


We rewrite our system as 
where 


Yt = Y{W -{N -2)-Y) for all t G M, 

Zi = Z{N +pY -Z) for all t G M, 

Wt = W{qZ -qN + q + N -W) for alH G 

Ct = 5(C) 


(7.4) 


(7.5) 



and g{C) = 


Y{W -{N -2)-Y) 
Z{N + pY -Z) 
W{qZ -qN + q + N -W) 


Since the system (17.5p is cooperative, the following comparison principle holds: 

Y{t) 


Lemma 7.1. Let C(t) = 
such that 



and C(t) = 


Z{t) 

W{t) 


be two nonnegative solutions of (f7.5D 


Y{to)>Y{to), Z{to)>Z{to), W{to)>W{to) 


for some to gM.. Then 


Y{t) > Y{t), Z{t) > Z{t), W{t) > W{t) for all t > to. 


From (17.11) and (17.21) we have Z > N and W > N. Therefore there are only two equilibria of 
(EH which satisfy Z > N and IT > A^, namely 


Ci= IV 
V N + q 


and (2 


2+9 


l-pq 

p{2+q) 

l-pq 

\ A^-2 + +T 

\ l-pq 


N + 


\ 
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Lemma 7.2. C ,2 is asymptotically stable. 
Proof. The linearized matrix at (2 is 


M 


-Y 2 0 Y 2 

pZ2 —Z2 0 

0 qW2 —W2 


The characteristic polynomial of M is 


P{X) = det(A/ — M) = A 3 + aA 2 + /3A + (1 — pq)'^ 


where 

a. = Y2 + Z 2 + W 2 

P = Y 2 Z 2 + ^ 2^2 + ^ 2^2 

7 = Y2Z2W2. 

Since a, /3, 7 > 0 and pq < 1, we have -P(A) > 0 for all A > 0. If -P has three real roots then they 
are all negative, so C 2 is asymptotically stable in this case. It remains to consider the situation 
where P has exactly one real root. Let Ai G M and A 2 , A 3 G C \ M be the roots of P. We claim that 
Re(A 2 ) = Re(A 3 ) < 0. We need to show P{—a) = —[3a + (1 — pq)j < 0, that is, fta > (1 — pq)^. 
By AM-GM inequality we find 

a > 3^/y^Z^ and [3 > 


which yields a/3 > (1 — pq)'y. Hence (2 is asymptotically stable. 
Lemma 7.3. For all t G M, we have 


0 < Y{t) < 


2 q 


N < Z{t) <N + 


1 ’ 

1 — pq 

p{2 + q) 


N + q< W{t) <N-2 + 


1-pq ' 

2 + q 


1 — pq 


Proof. Since u'(0) = 0 and v{0) = 0 we obtain limt_!._oo Y(t) = lim^-^o = 0. Also 


N — 1 

u"{r) H- u'{r) = u^(r) for all r > 0 . 

r 


Using L’Hopital’s rule, we deduce that 


v! ( 

lim- = lim u"(r) = u"(0) 

r^-o r r^O 


□ 


(7.6) 

(7.7) 

(7.8) 


(7.9) 


which combined with (j7.9p yields 


«"( 0 ) 


u^’(O) ^ u'(r) 

—= hm- 

N r —>-0 r 
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Therefore 


rv^(r) 

lim Z (t) = lim —= N. 
t->—oo r—^o u'[r) 

We claim that there exists tj ^ —oo such that 

' Y{tj) < Y2, 

< Z{tj) < Z 2 , (7.10) 

, W{tj) < W 2 . 


Because limt^_oo T(t) = 0 and limt^_oo Z(t) = N, it remains only to prove the last part of (jT.lOp . 
Assume by contradiction that this is not true. Thus W > W 2 in {— 00 , to) for some to € M. Then, 
by taking to small enough we have 

Wt = W{qZ — qN + q + N — W) <0 in (— 00 , to)- 

Hence, W is decreasing in the neighbourhood of —00 and so, there exists i = limi_)._cx 3 W (t). Again 
using L’Hopital’s rule we have 


i = lim W (t) = lim — 
t ->—00 r—fO v'[r) 


= lim 
r—>-0 


= lim 
r—>-0 


= lim 

t—)• —00 

i + q 


rqu'^~^{r)u" (r) + u'^{r) 
v"{r) 

rqu'^~^{r)v'P{r) + (1 + g — qN)u'^{r) 
u'i[r) — ^^^^v'{r) 
qZ{t) + 1 + q - qN 


1 - 


jV-l 

W{t) 


1 - 


Af-1 ■ 


This yields i = N q < W 2 which contradicts our assumption that W > W 2 in a neighbourhood 
of — 00 . This proves the last part of (IT.lOjl . We then apply the Comparison Lemma l7.II on all the 
intervals [fj,oo) for j > 1 to obtain the upper bound inequalities in Lemma 17.31 In the same way 
we obtain the lower bound inequalities. □ 


Let K = [[(^ 1 , (( 2 ]] C By Lemma 17?^ we have uj{C) Y K. Since C 2 is asymptotically stable, K 

has no circuits. Also, by (I2.17j) we have 

div g{C) = —W + {q — 2)Z + {p — 2)Y + N + 2 — qN + < 0 in A. 


Using Theorems 13.41 and 13.51 we deduce that uj{C) reduces to one of the equilibria Ci or ^ 2 - If 
C(t) — >■ Cl 9.S t —>■ 00 this implies in particular that Y{t) —0 as f —)• 00 . On the other hand, from 
the second equation of (|7.3I) we deduce T) > 0 in a neighbourhood of infinity which is impossible 
given that Y{t) >0 in M. Hence ({t) —>■ C 2 as t ^ 00 , that is 


hm X{t)=2 + P^^, 
t^oo 1 — pq 

2 + 0 

lim Y{t) = - -, 

t^Qo 1 — pq 

ii„,z(f) = jv+hl±h, 

t^oo 1 — pq 

lim W(t) = N -2 + 

t^oo ^ — pq 


26 























Using {X{t),Y{t), Z[t),W{t)) , we have 


lim 

r—^oo 


v{r) 

r 1-pg 


lim (y(t)Z'^(t)VU(i))p9-i 


{Y2ZIW2)^ 

{2Yq){N{l -pq) + p{2 + q))'^{N{1 - pq) + q{2p + 1))- ^ 


(1 - pq)'^+i 


And 


lim 


u{r) 

2+2p-pq 

T 1-pg 


{Y{t)Z‘i{t)W{t))^^ 

^4™ x{t)z{t) 

{Y2Z^W2)^ 

X 2 Z 2 

(2 + q){N{l -pq) + p{2 


ATf-\ 


2p+2—pg 

{l-pq) p {2 + 2p-pq) 


pg-1 

p 
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